home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / slar1v.z / slar1v
Encoding:
Text File  |  2002-10-03  |  5.5 KB  |  199 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAARRRR1111VVVV((((3333SSSS))))                                                          SSSSLLLLAAAARRRR1111VVVV((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLAR1V - compute the (scaled) r-th column of the inverse of the
  10.      sumbmatrix in rows B1 through BN of the tridiagonal matrix L D L^T -
  11.      sigma I
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SLAR1V( N, B1, BN, SIGMA, D, L, LD, LLD, GERSCH, Z, ZTZ,
  15.                         MINGMA, R, ISUPPZ, WORK )
  16.  
  17.          INTEGER        B1, BN, N, R
  18.  
  19.          REAL           MINGMA, SIGMA, ZTZ
  20.  
  21.          INTEGER        ISUPPZ( * )
  22.  
  23.          REAL           D( * ), GERSCH( * ), L( * ), LD( * ), LLD( * ), WORK(
  24.                         * ), Z( * )
  25.  
  26. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  27.      These routines are part of the SCSL Scientific Library and can be loaded
  28.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  29.      directs the linker to use the multi-processor version of the library.
  30.  
  31.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  32.      4 bytes (32 bits). Another version of SCSL is available in which integers
  33.      are 8 bytes (64 bits).  This version allows the user access to larger
  34.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  35.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  36.      only one of the two versions; 4-byte integer and 8-byte integer library
  37.      calls cannot be mixed.
  38.  
  39. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  40.      SLAR1V computes the (scaled) r-th column of the inverse of the sumbmatrix
  41.      in rows B1 through BN of the tridiagonal matrix L D L^T - sigma I. The
  42.      following steps accomplish this computation : (a) Stationary qd
  43.      transform,  L D L^T - sigma I = L(+) D(+) L(+)^T, (b) Progressive qd
  44.      transform, L D L^T - sigma I = U(-) D(-) U(-)^T, (c) Computation of the
  45.      diagonal elements of the inverse of
  46.          L D L^T - sigma I by combining the above transforms, and choosing
  47.          r as the index where the diagonal of the inverse is (one of the)
  48.          largest in magnitude.
  49.      (d) Computation of the (scaled) r-th column of the inverse using the
  50.          twisted factorization obtained by combining the top part of the
  51.          the stationary and the bottom part of the progressive transform.
  52.  
  53.  
  54. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  55.      N        (input) INTEGER
  56.               The order of the matrix L D L^T.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAARRRR1111VVVV((((3333SSSS))))                                                          SSSSLLLLAAAARRRR1111VVVV((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      B1       (input) INTEGER
  75.               First index of the submatrix of L D L^T.
  76.  
  77.      BN       (input) INTEGER
  78.               Last index of the submatrix of L D L^T.
  79.  
  80.      SIGMA    (input) REAL
  81.               The shift. Initially, when R = 0, SIGMA should be a good
  82.               approximation to an eigenvalue of L D L^T.
  83.  
  84.      L        (input) REAL array, dimension (N-1)
  85.               The (n-1) subdiagonal elements of the unit bidiagonal matrix L,
  86.               in elements 1 to N-1.
  87.  
  88.      D        (input) REAL array, dimension (N)
  89.               The n diagonal elements of the diagonal matrix D.
  90.  
  91.      LD       (input) REAL array, dimension (N-1)
  92.               The n-1 elements L(i)*D(i).
  93.  
  94.      LLD      (input) REAL array, dimension (N-1)
  95.               The n-1 elements L(i)*L(i)*D(i).
  96.  
  97.      GERSCH   (input) REAL array, dimension (2*N)
  98.               The n Gerschgorin intervals. These are used to restrict the
  99.               initial search for R, when R is input as 0.
  100.  
  101.      Z        (output) REAL array, dimension (N)
  102.               The (scaled) r-th column of the inverse. Z(R) is returned to be
  103.               1.
  104.  
  105.      ZTZ      (output) REAL
  106.               The square of the norm of Z.
  107.  
  108.      MINGMA   (output) REAL
  109.               The reciprocal of the largest (in magnitude) diagonal element of
  110.               the inverse of L D L^T - sigma I.
  111.  
  112.      R        (input/output) INTEGER
  113.               Initially, R should be input to be 0 and is then output as the
  114.               index where the diagonal element of the inverse is largest in
  115.               magnitude. In later iterations, this same value of R should be
  116.               input.
  117.  
  118.      ISUPPZ   (output) INTEGER array, dimension (2)
  119.               The support of the vector in Z, i.e., the vector Z is nonzero
  120.               only in elements ISUPPZ(1) through ISUPPZ( 2 ).
  121.  
  122.      WORK     (workspace) REAL array, dimension (4*N)
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136. SSSSLLLLAAAARRRR1111VVVV((((3333SSSS))))                                                          SSSSLLLLAAAARRRR1111VVVV((((3333SSSS))))
  137.  
  138.  
  139.  
  140. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  141.      Based on contributions by
  142.         Inderjit Dhillon, IBM Almaden, USA
  143.         Osni Marques, LBNL/NERSC, USA
  144.  
  145.  
  146. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  147.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  148.  
  149.      This man page is available only online.
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  
  192.  
  193.  
  194.  
  195.                                                                         PPPPaaaaggggeeee 3333
  196.  
  197.  
  198.  
  199.